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FOREWORD 


This fourth volume of Cosmic-Ray Results is a com- 
pilation of data for the past nine years, assembled by the 
Department of Terrestrial Magnetism, Carnegie Institution 
of Washington. The present volume contains data from 
Carnegie Institution of Washington’s Compton-Bennett 
cosmic-ray ionization chambers at Fredericksburg, Virginia; 
Huancayo, Peru; Christchurch, New Zealand; and God- 
havn, Greenland. These stations have been operated and 
maintained with care by many local personnel who have 
an abiding faith that long-term observations of cosmic 
phenomena are useful. 

Men are short-lived on the cosmic scale—yielding their 
personal abundance and contributions over thirty years 
or so in any selected area of their interest. Can a cos- 
mologist add significantly to understanding during his life- 
time of only one part in 3 X 10° of the cosmic cycle, which 
is now thought to be about 80 billion years? The answer is 
clearly yes, and emphatically so. Not all the meaning or 
physical understanding certainly. But parts of it, yes. 
Enough parts, indeed, to augment our understanding of 
the nature of the universe about us, and to map some ways 
to continue the search for knowledge. 

From the data in the present volume and those pre- 
viously compiled, an apparent systematic variation with a 
period of two sunspot cycles has been discovered in the 
diurnal variation. This variation is in phase with the cycli- 
cal reversal of the magnetic field of our nearest star, the 
Sun.* A longer time of observation—more than one and 
a half cycles—is desirable to reaffirm this effect which, 
though not understood, is doubtless related to the Sun’s 
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general magnetic field and that in the spiral streams from 
the Sun. 

Volume XXII, mostly numbers on approximately 300 
pages, will add to the prior volumes. Other generations 
will have them for their consideration; hence, we are proud 
to submit them to the indelible page. They are irreplace- 
able and irrevocable history; records, fallible in technical 
acumen as viewed by future scientists, perhaps, but none- 
theless the best we know how to obtain within our means, 
intellectual and financial. The metaphor is not strained in 
comparing these less dramatic records to Darwin’s careful 
observations, and Captain Cook’s log of exploration—the 
former to give rise to the idea of organic evolution, the 
latter to yield the conviction of the relation of the con- 
tinents one to the other and of unsuspected islands between. 

This volume is for the specialist. We make no apologies. 
We have the faith that among the community of specialists 
there will be one and perhaps many who will proceed 
deeply and learn how to dig significant signals out of what 
is apparent noise and to teach, if not us, new generations 
new qualities of the universe we strive to understand and 
appreciate, and turn their learning toward the benefit of 
mankind. Few, outside this Department, know the metic- 
ulous care which Liselotte Beach and Scott Forbush exer- 
cised in preparing this volume. For two decades and until 
very recently, Mrs. Beach carried out the compilation en- 
tirely “by hand.” Scott Forbush has lent a lifetime of wis- 
dom to these volumes, of which Volume XXII will be, we 
trust, not the last. 


ELLis T. BOLTON, Director 


PREFACE 


The Compton-Bennett ionization chambers which pro- 
vided data for this volume were operated and maintained 
with the unselfish enthusiastic cooperation of several organ- 
izations and individuals. For making available the records 
from equipment belonging to the Huancayo Magnetic 
Observatory we are grateful to Alberto A. Giesecke, Jr., 
Director, and to the staff of the Instituto Geofisico del 
Peru. 

We gratefully acknowledge our debt to the Director 
of the U. S. Coast and Geodetic Survey, Environmental 
Science Services Administration, for providing housing for 
the cosmic-ray meter at its Fredericksburg Geomagnetic 
Center. We express our appreciation to the Chief and the 
staff of the Fredericksburg Geomagnetic Center for their 
efficient operation and maintenance of the equipment. We 
acknowledge the cooperation of the Director, Magnetic 
Survey Division, Department of Scientific and Industrial 
Research of New Zealand, for making available the facil- 
ities of the Geophysical Observatory at Christchurch; and 
to the Director and staff of that observatory we are par- 


ticularly grateful, not only for efficient operation and 
maintenance of the meter, but also for scaling all the 
hourly cosmic-ray records for that station from April 1936 
to December 1961. We thank the Director of Det Danske 
Meteorologiske Institut, and the Director and staff of its 
Geofysisk Observatorium at Godhavn, Greenland, for the 
effective operation and maintenance of the meter there. 
Grateful acknowledgment is made to the Government of 
Denmark for transportation of supplies to Godhavn. 

On July 1, 1959, the cosmic-ray meters at Christchurch, 
New Zealand, and at Godhavn, Greenland, were turned 
over to, respectively, the Department of Scientific and 
Industrial Research and the Danske Meteorologiske In- 
stitut. 


LISELOTTE BEACH 

ScoTtr E. ForBUSH 

Department of Terrestrial Magnetism 
March 1, 1969 
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Cosmic Ray Results from Huancayo, Peru; 


Fredericksburg, Virginia; Christchurch, New Zealand; 
and Godhavn, Greenland 


INTRODUCTION 


Volume XIV and volume XX of Researches of the 
Department of Terrestrial Magnetism (Carnegie Institution 
of Washington Publication 175) contained cosmic-ray data 
from Huancayo, Cheltenham, Christchurch, and Godhayn. 
Volume XIV contained bihourly mean values for Huancayo 
and summaries of data for all four stations through 1946. 
Volume XX extended the summaries from Huancayo, 
Cheltenham, and Christchurch through 1955 and from 
Godhavn through 1950, and in addition included not only 
bihourly mean values from Huancayo for the period Jan- 
uary 1946 through December 1955 but also bihourly mean 
values from Cheitenham for the period March 1937 
through December 1955. 

Volume XIV described a typical installation, the instru- 
ment used, the absolute values of ionization, and proce- 
dures for reduction of the data. It also gave the geographi- 
cal coordinates, the altitude, and the geomagnetic latitude 
of the stations. The Fredericksburg Magnetic Observatory 
is located at geographic latitude 38°12’N and longitude 
77°22'W at an elevation of 69 meters above sea level. 


Volume XXI included data subsequently obtained 
through 1959 from Huancayo and Cheltenham-Fredericks- 
burg. The station at Cheltenham was moved to Fredericks- 
burg on October 5, 1956; care was taken to ensure that 
there is no discontinuity in the published values at the 
time of this transfer. Bihourly means and summaries for 
Huancayo and Cheltenham-Fredericksburg were included. 
Summaries for Christchurch were for the period January 
1956 through June 195°. Bihourly means and summaries 
for Godhavn are for the period January 1954 through July 
1959. Also included are bihourly means and summaries 
from Universidad Nacional de México, Instituto de 
Geofisica, for the IGY period July 1957 through Decem- 
ber 1958. 

The present volume contains bihourly means and sum- 
maries: for Huancayo and for Fredericksburg, January 
1960 to December 1968, and for Godhavn, January 1951 
to December 1953, previously unpublished; for Christ- 
church only summaries, January 1959 to December 1961, 
are included. 


REDUCTION OF DATA 


As indicated in the Introduction, volumes XIV and XX 
described some details of the procedures for reducing the 
data at each of the stations, with further details in volume 
XXI. These include a table of corrections to the 24-hour 
and 12-hour harmonic coefficients for the pressure-cor- 
rected ionization at Huancayo to compensate for the bi- 
hourly mean values of pressure having been erroneously 
centered one-half hour late from January 1946 to Decem- 
ber 1958. 

In a paper entitled “Diurnal Variation in Cosmic-Ray 
Intensity, 1937-1959, at Cheltenham (Fredericksburg), 
Huancayo, and Christchurch” by Scott E. Forbush and D. 
Venkatesan, tables are given for the annual mean 24-hour 
and 12-hour waves at these three stations. In these tables, 
the values of these harmonic coefficients for Huancayo for 
the period January 1946 through December 1958 are those 
that would be obtained from properly applied pressure 
corrections during that interval. This paper was published 
in the Journal of Geophysical Research, volume 65, pages 
2213-2226, August 1960. It also contains a critical ex- 
amination of the errors that can arise from inevitable fric- 
tion (and imperfect temperature compensation) in record- 
ing barographs and microbarographs. 


Results of a more recent investigation of errors due to 
friction in the barograph in the annual mean harmonic co- 
efficients for barometric pressure at Huancayo are de- 
scribed in a memorandum at the end of this section under 
the heading: Barometric pressure corrections to the 
cosmic-ray ionization at Huancayo and the change in ion- 
chamber sensitivity after about 1955. This memorandum 
contains tables of the annual mean values of the harmonic 
coefficients for ionization at Huancayo, with and without 
the pressure corrections derived from the investigation, 
together with corrections for a change in sensitivity for the 
ionization after about 1955. 

At the several stations the argon pressure in the ioniza- 
tion chamber has remained constant through 1968 at the 
values given in volume XXI. 

It is important to note that all data in the present volume, 
as well as in volumes XIV, XX, and XXI, are corrected 
only for barometric pressure. No other meteorological 
corrections have been applied. Consequently, seasonal 
waves are present in the data for all stations except perhaps 
Huancayo. In addition, the data at all stations except 
Cheltenham indicate slow changes in the balance current. 
These instrumental drifts and seasonal variations are de- 
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scribed in “World-Wide Cosmic-Ray Variations, 1937- 
1952,” Journal of Geophysical Research, volume 59, pages 
525-542, December 1954. 

It now appears that the rate of this instrumental drift 
(changing balance current) at Huancayo changes after 
adjustments of the uranium rod are made to alter the com- 
pensation, or balance, for the cosmic-ray ionization. 

At Huancayo this compensation was changed by ad- 
justing the uranium rod on the following dates: March 22, 
1954, October 2, 1957, November 10, 1964, and Decem- 
ber 28, 1966. As described in earlier volumes certain con- 
stant values were arbitrarily added (in the scaling proce- 
dure) to the actual departures from balance to ensure that 
negative tabular values would seldom occur. These addi- 
tive values in units of 0.1 percent are as follows: 


+ 22 from June 1936 to March 1937 

+ 45 from March 1937 to February 1941 

+ 60 from March 1941 to March 22, 1954 
+125 from March 22, 1954, to October 2, 1957 
+100 from October 2, 1957, to December 1968 


However, it should be noted that on all tables for the 
period after February 1941 the headings indicate that +60 
(in units of 0.1 percent) has been added to the tabulated 
values. The apparent discrepancy between this value of 


+60 and the value +125 used from March 22, 1954, to 
October 2, 1957, as shown above, arises from the fact that 
the value +125 includes an additive value of +65 to com- 
pensate the change of —65 (in units of 0.1 percent) due 
to adjusting the uranium rod on March 22, 1954. Thus, 
for the period March 22, 1954, to October 2, 1957, the 
tabulated values are the estimated departures from balance 
(with +60 added) that would have prevailed (as nearly 
as can be determined) had the uranium rod not been ad- 
justed on March 22, 1954. Similarly the actual change in 
the departure from balance due to the adjustment of the 
uranium rod on October 2, 1957, was —25 (in units of 0.1 
percent). Thus the value of +100 actually added (in the 
scaling process) to the departures’ from balance from 
October 2, 1957, to December 1968 results in estimated 
departures from balance (with +60 added) that would 
have prevailed, until November 10, 1964, had the uranium 
rod not been adjusted on March 22, 1954, and on October 
2, 1957. However, from November 10, 1964, to December 
1968 the value of +100 in units of 0.1 percent added to 
the actual departures from balance does not correct for 
the change due to adjusting the uranium rod on November 
10, 1964, and on December 28, 1966. Thus, after Novem- 
ber 10, 1964, the value +100 in units of 0.1 percent was 
added to the actual departures from balance and not the 
value +60 listed in the tables. 


CORRECTIONS TO FOURIER COEFFICIENTS FOR THE DIURNAL VARIATION AT 
CHRISTCHURCH, 1956 TO 1961, TO COMPENSATE FOR AN ERROR IN THE 


CORRECTIONS 


In October 1969 it was found that for Christchurch, 
from March 6, 1956, to December 31, 1961, the bi- 
hourly means of the barometric pressure corrections 
were centered one-half hour earlier than the uncor- 
rected bihourly means of the departures from balance 
to which the corrections were applied. 

The following Table A gives the consequent correc- 
tions Aa,, Ab, and Aas, Abs in units of 0.01%, to be 
added to the Fourier coefficients a,, 5; and ao, b» re- 
ferred to GMT and to 172.5° EMT, for the yearly mean 
24-hour and 12-hour waves derived from the published 
diurnal variations from 1956 through 1961: 


Table A. Corrections to Fourier Coefficients in 
Units of 0.01% 


Time origin Aa, Ab, Aa» Abs 

GMT —().73 +0.34 —1.65 —1.30 
172.5° EMT +0.68 —0.43 —1.92 —().82 
For b, —O and b. — 0 the times of maximum are at 


the hour zero. As a consequence of the small variation 


FOR PRESSURE 


among yearly mean values of the Fourier coefficients 
for the annual mean diurnal variation in barometric 
pressure shown in the following Table B for six arbi- 
trarily chosen years, the above corrections (derived 
from the six-year means at the bottom of Table B) are 
essentially constant for the years 1956 to 1961: 


Table B. Fourier Coefficients, Origin GMT, in Units of 0.01 
Millibar for Each of Six Years in the Annual Mean Diurnal 
Variation of Barometric Pressure at Christchurch 


Year a by a» bs 
1954 —22 —34 +31 —52 
1955 —21 —47 +28 —48 
1956 —16 —24 +28 —47 
1960 —22 —50 +30 =—55 
1961 —26 —41 +33 —55 
From May 

1965 —24 —42 +29 ==) 1 
Mean 21.7 —39.8 +29.9 —51.4 


The barometric coefficient 6 —= 0.135% /millibar was 
used to obtain the values in Table A from the means 
in Table B. 


BAROMETRIC PRESSURE CORRECTIONS TO THE COSMIC-RAY IONIZATION AT 
HUANCAYO AND THE CHANGE IN ION-CHAMBER SENSITIVITY AFTER ABOUT 1955 


ScoTT E. FORBUSH 


From June 1936 through February 1958 all (published 
and otherwise available) bihourly values of ionization at 
Huancayo have been corrected using pressure, P,, from the 
ordinary Richard barograph and barometric coefficient 
Bo = —0.30%/mm Hg. This value of 8B, = —0.30%/mm 
Hg, for the ordinary barograph, was obtained from data 
during the period June to September 1937. For Huancayo, 
where daily means of pressure vary only a few mm Hg 
during any month, a special procedure used to determine 
B has been described!?. In this procedure, intervals of 
about 4 consecutive days are selected during which the 
daily means of pressure increase (or decrease) monotoni- 
cally over a minimum range of at least 1.5 or 2.0 mm Hg. 
For a number (the more the better) of such intervals of 
increasing (or decreasing) pressure the mean departures, 
from the average pressure for all days and all intervals, are 
obtained for interval days 1, 2, 3 and 4. Corresponding 
departures for ionization uncorrected for pressure are also 
derived, and plotted as function of the pressure departures 
to determine the barometric coefficient, 8. This procedure 
provides, within sampling uncertainties, the correct value 
of 8 if the daily means of pressure are not subject to appre- 
ciable statistical uncertainty. Using this procedure, and 
data for the period January 1954 to December 1955, the 
value 8, = —0.306%/mm Hg was obtained’ for the 
ordinary barograph pressure. 

Beginning March 1958 published bihourly values of 
ionization at Huancayo were corrected using pressure, P;,, 
from a microbarograph and _ the pressure coefficient, 
B = —0.30%/mm Hg, that had been obtained for the 
ordinary barograph. 

Table A lists the yearly mean harmonic coefficients, on 
75° WMT, for barometric pressure and for ionization, 
1(P), uncorrected for pressure. It should be noted that the 
coefficients for pressure in mm Hg are obtained by dividing 
the tabulated values by 30. From 1955 to 1966 two sets 
of coefficients are listed for ionization, 7(P), uncorrected 
for pressure. Those in parentheses have been corrected for 
an increase of 10%, relative to that before 1955, in ion- 
chamber sensitivity which is approximately that determined 
by methods to be described. 

Figure 1 shows that the 75° WMT times of maximum 
for the yearly 24-hour wave are correlated with those for 
the 12-hour wave from the ordinary barograph, 1937-1957, 
and from the microbarograph, 1958-1966. The slope, 1.90, 
of the line in figure 1 indicates that the variations in phase 
angle for the 24-hour wave are 0.95 times those for the 
12-hour wave or that the variations in the two phase angles 
are nearly the same. This indicates that these variations 
are likely due to variations from year to year in the fric- 
tion of the pen-recording ordinary barograph and in the 
pen-recording microbarograph. Since the times of maxi- 
mum for the 1937-1941 average 12-hour and 24-hour 
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Fig. 1. Correlation between times of maximum for 12-hour 


and 24-hour waves in pressure for different barographs at 
Huancayo. Tmax on 75° WMT. 


wave from the ordinary barograph were both earlier than 
for averages for later periods (see, for example, rows | and 
2 of table B) the coefficients in row 1 of table B were 
adopted for each year 1937 to 1966. Also, table B (rows 
1 and 8) shows that the time of maximum of the average 
(1937-1941) 12-hour pressure wave was 6"08™ later than 
the average (1937-1941) 12-hour wave in ionization, 
I(P), uncorrected for pressure. The phases of these 12- 
hour waves are thus almost exactly 180° apart, indicating 
that the 12-hour wave in uncorrected ionization is prob- 
ably due entirely to the 12-hour wave in pressure. Thus, 
for the averages, 1937-1941, C.[I(P)]/C.(P,) or 
bo [I(P)]/b.(P) were used to determine the pressure co- 
efficient, 8, = —0.33%/mm Hg, for the ordinary bar- 
ograph. This value of 8, = —0.33%/mm Hg and the 
Fourier coefficients, for pressure, in row | of table B give: 
Aa, = +8 and Ab, = +34 in units of 0.01%, for the 
respective corrections to be added to the values of 
a,U(P)] and 6,[I(P)] for each year, as indicated at the 
bottom of table C. These corrections were used to obtain 


Barometric Pressure P 
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Table A. Coordinates in 24-hour and 12-hour Harmonic Dials on 75° WMT for Yearly 
Mean Waves in Pressure and in Ionization Uncorrected for Pressure at Huancayo 


C—R Uncorrected for P 


u = 0.01% (1/30 mm Hg) u = 0.01% 
Year ay by de bo ay by ay bo 
1937 +7 +31 +11 —29 —19 —22 — 9 +31 
1938 +6 +32 +13 —28 —21 —24 —11 +30 
1939 +8 +29 +12 —28 —25 —27 —12 +31 
1940 +8 +30 +11 —29 —25 —29 —11 +33 
1941 +8 +32 +13 —29 —24 —30 —12 +33 
1942 +8 +32 +14 —29 —23 —32 —14 +32 
1943 +6 +30 +13 —28 <=25 —31 —13 +32 
1944 +5 +32 +15 —27 —15 —27 —13 +35 
1945 +6 +32 +14 —27 —19 —28 —13 +33 
1946 +6 +26 +10 —26 —24 —24 —10 +31 
1947 +8 +28 +10 —27 =25 —30 —14 +33 
1948 +6 +28 +13 —26 —20 —24 —13 +32 
1949 +7 +29 +11 —27 —22 —21 —10 +29 
1950 +6 +30 +11 —28 —24 =21 —12 +30 
1951 +7 +31 +12 —28 —27 —18 = 9 +32 
1952 +7 +29 +12 —27 —20 —16 —10 +31 
1953 +7 +31 +13 —27 =15 —13 —13 +32 
1954 +4 +29 +15 —25 — 4 —19 =I/5 +29 
1955 +3 +29 +16 —24 —12(—13)* —14(—-15) —10(—11) +31(+434) 
1956 +5 +29 +14 —25 —21(-—23) —21(—23) —11(—12) +30(+33) 
1957 +5 +29 +14 —25 —22(-—25) —23(—26)  —13(—14) +30(+33) 
1958 +4 +30 +14 —26 —22(—24)  —32(—35) —11(—12) +32(+35) 
1959 +4 +28 +15 —24 —24(—27) —31(-—34)  —12(—13) +33(+37) 
1960 +3 +28 +14 —22 —22(—24) -—32(-—36)  —11(-—12) +32(+36) 
1961 +1 +26 +14 —20 —15(—17) —25(—28)  —13(—14) +30(+33) 
1962 +7 +31 +17 —18 —24(—27)  —32(—36) —13(—14) +32(+35) 
1963 —1 +27 +15 —18 —17(—19)  —27(-—30) —13(—14) +31(+34) 
1964 +2 +28 +15 —21 —16(—18) —28(—31) —13(—14) +34(+38) 
1965 +4 +29 +13 —22 —14(—15)  —27(—30) —13(—14) +35(+39) 
1966 +5 +31 +13 —22 —22(—25) —32(—35)  —13(—14) +34(+38) 
Sum +162 +986 +397 —762 —608 —760 —360 +953 
Mean + 5.4 +29.5 +13.2 —25.4 —20.3 =25.3 —12.0 +31.8 
*Values in parentheses are without corrections for 10% increase of sensitivity. 
Note: 1937—Feb. 1958 from ordinary Richard barograph. Feb. 1958-1966 from microbarograph. 
Table B. Fourier Coefficients, at Huancayo, for 12- and 24-hour Waves: (A) in Barometric 
Pressure and (B) in Cosmic-Ray Ionization Uncorrected for Pressure, /(P). 
Tmax = 75° WMT of Maximum. 
(A) 
ay by C1 Tmax a2 bo cog Tmax 
Barograph Years mm Hg mm Hg mm Hg hrs min mm Hg mm Hg mm Hg _hrs min 
1. Ordinary 1937-1941 +0.247  +1.027 1.056 5 06 +0.400 —0.953 1.034 9 46 
2. Ordinary 1953-1957 +0.160  +0.980 0.993 5 23 +0.480 —0.840 0.967 9°59 
3. Simpson 1953-1957 +0.268 +0.944 0.981 4 57 +0.336 —0.933 0.992 9 40 
4. Microbarograph 1958-1966 +0.107 +0.957 0.963 5 34 +0.480 —0.713 0.860 10 08 
5. Simpson 1958-1966 +0.231 +0.936 0.964 5 05 +0.367 —0.930 1.000 9 43 
6. Simpson 1953-1966 +0.245 +0.939 0.970 5 01 +0.356 —0.931 0.997 9 42 
(B) © % % % 
7. I(P)* 1953-1966 —0.124 +0.318 0.341 3 43 
8. I(P) 1937-1941 —0.110 +0.316 0.335 338 
*Includes correction for 10% increase in sensitivity, 1955-1966. 
Note: Using cg from 8 and 1, 6 = —0.324%/mm Hg, for ordinary barograph. Using co from 7 and 6, Bs = —0.342%/mm 
Hg, for Simpson’s barograph. Using bs from 7 and 6, Bs = —0.342%/mm Hg, for Simpson’s barograph. Using cz from 7 and 6, 


Bs = —0.348%/mm Hg, for Simpson’s barograph. 
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the adopted yearly values of a, and 5, in table C for ioniza- 
tion at Huancayo, corrected for pressure. 

After these pressure corrections were adopted and used, 
data from Simpson’s barograph became available. Table 
B shows that the average, 1953-1966, 12-hour wave for 
Simpson’s barograph is exactly opposite in phase to that 
for uncorrected ionization. Table B also shows that the 
average, 1953-1957, 12-hour wave is practically the same 
as that for 1958-1966, and that the time of maximum for 
the first of these is 19 minutes earlier than the average 
12-hour wave from the ordinary barograph, while the sec- 
ond is 25 minutes earlier than that from the microbar- 
ograph. Table B also shows comparable differences in the 
times of maximum for the 24-hour waves. For the period 
1953-1966, a plot of Ty,,, for the 12-hour and for the 
24-hour waves, from Simpson’s barograph, for individual 
years showed no correlation. Also, all except two of the 
14 values of 7\,,, for the 12-hour wave were within +6 
minutes of the average T,,,, for the 14 years. Some of this 
variation is doubtless due to the fact that the yearly mean 
Fourier coefficients which were furnished were obtained 


from hourly means averaged for all days of the year. This 
of course is correct if data for the same number of days 
are available in each month. There is a considerable vari- 
ation in phase and amplitude of the monthly mean 12- and 
24-hour waves during the year at Huancayo. Thus, if the 
yearly mean waves are derived as first described, and data 
are missing or quite incomplete for some months (as was 
the case in some years), the annual mean waves will be 
in error. Table B shows for Simpson’s barometric pres- 
sure, Py, that By = —0.34%/mm Hg whether obtained 
from as(I[p)/ao(Pg), bo(Ip)/bo(Ps), or co(Ip)/cs(Pyg). 
The difference between this value By = —0.34%/mm Hg 
and the value By = —0.35% /mm Hg obtained from figure 
7 (discussed later) is hardly significant statistically. Thus 
there is no evidence for friction in Simpson’s barograph. 
Table D lists the Fourier coefficients for the yearly average 
24-hour and 12-hour waves and the times of maximum, 
from Simpson’s barograph at Huancayo. 

Figure 2 indicates that annual means of the 12-hour 
wave in uncorrected ionization were essentially constant 
from 1937 to about 1955, and also from about 1956 to 


Table C. Harmonic Coefficients for Pressure-Corrected Cosmic-Ray Ionization, 
Origin Standard Mean Time, in Units of 0.01% 


Huancayo* Fredericksburg} Christchurch Godhavn 
Year ay by a by ay by ay by 
1937 —i1 +12 —12 — 3 —11 — 3 
1938 —13 +10 —12 — 4 —15 — 4 
1939 —17 + 7 —14 —11 —14 —11 — 5 — 6 
1940 —17 + 5 —10 —11 —16 —14 — 5 6 
1941 —16 + 4 —10 —10 —14 —13 — 4 — 6 
1942 —15 + 2 — 8 —l1 — 9 —11 — 5 — 6 
1943 —17 + 3 —12 —12 —13 —11 — 6 — 9 
1944 — 7 + 7 — 7 —2 — 6 — 5 — 4 + 1 
1945 —11 + 6 — 8 — 7 — 8 —9 — § — 5 
1946 —16 +10 —10 —10 —13 — 8 — 6 —14 
1947 —17 + 4 — 5 —16 —10 —12 7 — 7 
1948 —12 +10 —14 —11 —12 — 5 — 8 —16 
1949 —14 +13 —14 — 9 —15 — 1 — 6 —17 
1950 —16 +13 —12 0 —14 —2 — 7 —11 
1951 —19 +16 —15 — 2 —16 0 — 9 — § 
1952 —12 +18 —15 0 —13 +2 — § 0 
1953 — 7 +21] —11 + 8 — 4 + 7 — 7 0 
1954 + 4 +15 0 + 6 — | + 5 0 4 
1955¢ — 4 +20 — 8 +4 — 7 + 8 — 5 + 2 
1956 —13 +13 —17 — 3 —13 + J — 9 — 4 
1957 —14 +11 —13 — 8 —15 — 1 — 9 — 4 
1958 —14 2 —10 — 8 —10 — 7 — 4 — 6 
1959 —16 ape — 9 —10 —10 7 
1960 —14 + 2 —10 —10 —10 6 
1961 — 7 + 9 — 7 — 3 — 6 6 
1962 —16 = a —12 —13 
1963 — 9 + 7 —11 — 7 
1964 — 8 + 6 —10 — 9 
1965 — 6 arid — 4 — 2 
1966 —14 + 2 — 9 — 7 
Sum — 368 +260 — 309 —181 — 275 —113 —116 —115 
Mean —12.3 +:8:7 —10.3 — 6.0 —11.0 — 4.5 — 5.8 — 5.8 
No. of ; 
years 30 30 25 20 


“For Huancayo the pressure corrections +8 for a; and +34 for b; (in units of 0.01% ). were used for each year. 
+Cheltenham, Jan. 1937—Sept. 1956; Fredericksburg, Oct. 1956-1966. 
_ }For Huancayo, 1955-1966, the values of a, and b, for ionization uncorrected for pressure were corrected for the 10% 
increase of sensitivity (table A) before applying the pressure corrections. 
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Table D. Yearly Mean Harmonic Coefficients in mm Hg from Simpson’s Barograph at Huancayo 


Dinax TL sivax 
75° WMT 75° WMT 

Year ay by hm ae be hm 
1953 0.254 0.932 4 0.301 —0.914 9 36 
1954 0.286 0.913 4 0.248 —0.937 9 30 
1955 0.220 0.946 5 08 0.369 —0.926 8 43 
1956 0.271 0.971 4 58 0.332 —0.967 9 38 
1957 0.311 0.957 4 48 0.432 —0.919 9 50 
1958 0.208 1.011 5 0.385 —0.987 9 43 
1959 0.257 0.886 4 0.363 —0.933 9 43 
1960 0.183 0.895 5 0.305 —0.853 9 39 
1961 0.217 0.901 5 0.367 —0.902 9 44 
1962 0.234 0.957 5 0.385 —0.958 9 44 
1963 0.224 0.894 5 04 0.342 —0.915 9 41 
1964 0.225 0.935 5 06 0.378 —0.932 9 44 
1965 0.232 0.967 5 06 0.363 —0.941 9 42 
1966 0.303 0.979 4 0.417 —0.947 9 48 
Sum 3.425 13.144 4.987 — 13.031 
Mean 0.245 0.939 5 02 0.356 — 0.931 9 42 


1966. From about 1955 to 1966, in figure 2, the amplitude 
of the average 12-hour wave was about 11% greater than 
from 1937 to about 1954. 

The 12-hour harmonic dial in figure 3 shows the waves 
for ionization uncorrected for pressure for each year, 
1937-1966, the average for 1937-1954, and for 1955- 
1966, and the so-called probable error circles for single 
years and for the mean of the two groups of years. From 
the latter error circles the probability that the difference 
between the means of the two groups arises from sampling 
errors is about 10° so the difference is statistically sig- 
nificant. Figure 3 shows that the time of maximum of the 
1937-1954 average 12-hour wave is practically identical 
with that for the 1955-1966 average. The time of maxi- 
mum for both these average 12-hour waves in ionization 
uncorrected for pressure is almost exactly 6 hours earlier 
than that for the 12-hour wave in barometric pressure from 
Simpson’s barograph. The amplitude of the average vector, 
*{955-1966, in figure 3, is about 12% (essentially the same 
as in figure 2) greater than the average for 1937-1954. 
This accounts for the result in figure 4 which shows that 
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Fig. 2. Annual means, 1937-1966, at Huancayo (A) for the 
component, 1, perpendicular and (B) for that parallel, P, to 
the average, 1937 to 1966, vector in the 12-hour harmonic dial 
for ionization uncorrected for pressure. (C) indicates cumu- 
lated annual values of [P — 35.0]. 
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Fig. 3. Harmonic dial on 75° WMT for yearly means of 12- 
hour wave in cosmic-ray ionization uncorrected for pressure 
at Huancayo 1937 to 1966. 


the barometric coefficient, By, (using the 12-hour waves 
in pressure and uncorrected ionization) for Simpson’s 
barograph is about 12% greater after about 1956 than be- 
fore. For figure 4 the values of b:(/») in parentheses in 
table A were used after 1954 (i.e., before correction for 
the 10% increase in ion-chamber sensitivity). Since a 
similar increase is found later for the barometric coefficient 
based on variations of daily means of pressure and of 
ionization uncorrected for pressure (and uncorrected for 
ion-chamber sensitivity) it is certain that the change of By 
shown in figure 4 is due to a change in ion-chamber sensi- 
tivity and not to some real change in the 12-hour wave in 
cosmic-ray intensity. 

Figures 5 and 6 show that the ratio of changes in daily 
means from the microbarograph and from Simpson’s 
barograph may be taken as equal. In addition, the corre- 
lation coefficient, r — +0.990, in figure 5, shows that the 
statistical uncertainty in daily means from Simpson’s 
barograph and from the microbarograph is quite small. 

Figure 7 indicates, for selected intervals (with complete 
daily data) in the period 1953-1955, the value By = 
—0.35% /mm Hg derived from variations of daily means 
of pressure and of ionization uncorrected for pressure 
using the special procedure described earlier and in refer- 
ences | and 2. This value of By = —0.35%/mm Hg as 
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showing change in barometric coefficient for Simpson’s_ bar- 
ograph due to increase in ion-chamber sensitivity after about 
1956. b,(U/r) and b,(Ps) are the yearly mean sine coefficients 
for the 12-hour wave (on 75° WMT) in ionization uncor- 
rected for pressure and for pressure from Simpson’s barograph, 
respectively. 


mentioned earlier does not differ significantly (in the 
statistical sense) .from that, By = —0.34% /mm Hg shown 
in the footnote to table B. 

Figure 8 shows the value 8B), = —0.38%/mm Hg de- 


rived from daily means, for selected 4-day intervals, dur- 
ing the period November 1964 to March 1966, using pres- 
sure from the microbarograph. Figures 5 and 6 indicate 
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Fig. 5. Daily mean pressure, Py, from microbarograph as func- 
tion of that, Ps, from Simpson’s barograph for 50 days with 
complete data at Huancayo, January 1964—April 1965. 
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Simpson’s barograph at Huancayo, 1964-1965, 


that the value of Bs for the period November 1964 to 
March 1966 would be practically the same, i.e., Bg = 
—0.38%/mm Hg for Simpson’s barograph. The ratio, 
—0.38/—0.34, indicates that the ion-chamber sensitivity 
was about 12% greater during the period November 1964 
to March 1966 than during the period 1953-1955, and 
that the similar increase in sensitivity deduced earlier from 
the 12-hour waves was not due to a real increase in the 
12-hour wave in cosmic-ray ionization. 

Figure 9 shows the regression lines derived from the 31 
daily means of pressure, P,, from the ordinary barograph 
and from the microbarograph, P,,, for March 1958. Fig- 
ure 10 shows the regression lines for the changes in pres- 
sure derived from 7 selected intervals of 3-S days with 
pressure increasing (or decreasing). 

The slopes of the regression lines in figure 9 (or 10) 
determine the correlation coefficient, r(P,, Py) = +0.94, 
between P, and P,,. Figure 5 shows r(Px,, P3,) = +0.99 
for the correlation between daily means, Pg, from Simp- 
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Fig. 7. Pressure coefficients, 6, for ion chamber at Huancayo, 
from data for several 4-day intervals, 1953-1955, during each 
of which daily means of pressure from Simpson’s barograph 
were increasing (or decreasing). 
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Fig. 8. Pressure coefficient for ion chamber at Huancayo from 
data for 40 four-day intervals, November 1964 to March 1966, 
with pressure increasing (or decreasing) using microbarograph. 


son’s barograph and those, P,,, from the microbarograph. 
Since {1 — rf? (Pg, Py)} = 0.01 and {1 — P(P,, Par) } 
= 0.11, the variance of residuals for P;, computed from 
P, is 11 times that for P;, computed from Ps; consequently 
the error in daily means of Pj, (or Py) is practically 
negligible compared to that in daily means of P,. Thus in 
figure 9 (or 10) the “true” value of P” (in terms of P),) 
is essentially that given by: P” = P,,/0.84 whereas when 
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Fig. 9. (A) Group means which determine the regression of 
daily means, P., for pressure from microbarograph on those, 
P., from ordinary barograph, and (B) for regression of P. on 
P,,. Group means for largest P. (or Px) are from 7 daily 
averages: the other 5 group means are from 6 daily averages, 
Huancayo, March 1958. 
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Fig. 10. Pressure changes, 6Ps,, from microbarograph as func- 
tion of those, 6P., from ordinary barograph from data (March 
1958) for 4 intervals of 3-5 days with increasing pressure, 
combined (after change of sign) with data for 3 such intervals 
with decreasing pressure as determined by P. for (A) and by 
P, for (B). Coordinates in units of 0.1 mm Hg for Huancayo, 
March 1958. 


P, is (etroneously) taken free of error (i.e., all error in 
Pj,) an erroneous value, P”, computed from P), is given 
by P? = P,,/0.75; consequently P? = 1.12 P’. Thus 
in deriving the barometric coefficient, 8,, for the ordinary 
barograph when pressures, P,, from the ordinary barograph 
were assumed free of error, the resulting value of £,, say 
BY", relative to 8, (from the microbarograph) is given by 
BY = Bx/1.20 or from 8” the value By = 1.20 8%. The 
value B” = —0.306% /mm Hg was derived, as mentioned 
earlier, from P, using data for the period January 1954 to 
December 1955. Thus By, = —1.20 X 0.306 = 
—0.342% /mm Hg, which agrees with the value By in 
figure 4 (see also table B) for the period 1953-1956. This 


value By, = —0.342 is also in reasonable agreement with 
the value Bs = —0.354 shown in figure 7, and with the 
value B,, = —0.385 shown in figure 8, when this is re- 


duced 10% (to correct for change in ion-chamber sensitiv- 
ity) to give By, = —0.350%/mm Hg. 

Finally, it should be noted that the most reliable values 
for the pressure corrections Aa,, Ab,, Aa, Ab, to be added 
to each yearly mean of corresponding Fourier coefficients 
from ionization uncorrected for pressure are those ob- 
tained from the pressure coefficient By = —0.34% /mm Hg 
using pressure from Simpson’s friction-free barograph (see 
footnotes of table B) and the Fourier coefficients for Py 
given in row 6 of table B. Thus in units of 0.01% these 
additive corrections (using subscript s to indicate they are 
determined from Simpson’s barometric data) are: 


Aayg = 0.245 X 34 = +8 
Abig = 0.939 X 34 = +32 
Aasg = 0.356 X 34 = +12 
Abog = —0.931 X 34 = —32 


The values obtained (before Simpson’s barometric data 
were available) from the ordinary barograph, using B, = 
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—0.33%/mm Hg and the Fourier coefficients for pressure 
in row 1 of table B, as described earlier, and used to ob- 
tain a, and 5, in table C for Huancayo are, in units of 
0.01%: 


Adi, = +8 
Abo = +34 


For the 12-hour wave the corrections that would be 
similarly obtained from the ordinary barograph are, in 
units of 0.01%: 


Ag3 = 813 
Abs, — =31 


However, the values Aaj, and Abs, from Simpson’s bar- 
ograph are undoubtedly more reliable and should be used 
for correcting the 12-hour wave in ionization. 

Finally it may be noted that the variances, s*[a,(Ip)] 
and s?[b,(Jp)] for the values of a, and hb, at Huancayo 
uncorrected for pressure as given in table A (not using 
values in parentheses) are: 


or M? = s?(a,) + s2(bo) = 4.8 for the two-dimensional 
variance (M = Bartels’ expectancy). The corresponding 
values for the pressure corrected values, s?la,(/,)], and 
s*[b,(1,)], derived from the years 1937-1957 using for 
each year the Fourier coefficients obtained from data from 
the ordinary barograph and 8, = —0.33%/mm Hg, are: 


s?[ay 
s°Lb,(1,)] = 4.4 
M2 = s?[a,(I,) + s2b,(1,)] = 7.0 


Thus the two-dimensional variance, M?, is 50% greater 
after correcting for pressure than before correcting for 
pressure. It follows that M? = 4.8 if the same value of 
Aa, and of Abs is used for the pressure correction for 
each year. Thus for the 12-hour wave the use of constant 
corrections for pressure is justified, since M? = 7.0 when 
the actual corrections for each year (from the ordinary 
barograph) are used; figure | indicates that these constant 
corrections are certainly more reliable than the actual 
yearly corrections deduced from the microbarograph pres- 
sure. 

For the years 1937-1957 the variances of the yearly 
pressure corrections, Aa,, and Ab,,, in units of 0.01%, 


from the ordinary barograph using 8, = —0.33%/mm Hg, 
are: 


s2(Aa,,) = 2.35 
$(Ab,,) = 3.25 
M,2 = 5.60 


For the years 1953-1966 the variances of yearly cor- 
rections, Aa,, and Ab,, in units of 0.01%, from Simpson’s 
barograph using By = —0.34%/mm Hg are: 


s?(Aa,,) = 1.64 
s°(Ab,,) = 1.64 
Mg = 3.68 


Since there are certainly genuine variations in the yearly 
values for the 24-hour wave in cosmic-ray ionization, the 
above values of M? cannot be usefully compared with M? 
for the 24-hour yearly waves uncorrected for pressure. 
However, the values Aa,, = +8 and Ab,, = +34, men- 
tioned earlier, were used for Huancayo for each year. 
Deviations of yearly values from a 25-year mean have been 
derived for Cheltenham-Fredericksburg, Christchurch, and 
Huancayo. After these deviations were normalized (for 
phase and amplitude) to the deviations for Cheltenham- 
Fredericksburg, an analysis of variance gave for each sta- 
tion the residual variance 6.0 (with normalized deviations 
in units of 0.01%) for a single year. Since the amplitude 
normalization factor for Huancayo (to Cheltenham) was 
0.91, the residual variance for Huancayo is 7.3 (deviations 
in units of 0.01%), which may be compared with 6.0 for 
the residual variance for Cheltenham-Fredericksburg and 
Christchurch for which amplitude normalization factors 
were both unity. Thus it appears that the use of a con- 
stant pressure correction Aa, and Ab, for each year is also 
justified for Huancayo. 


REFERENCES 


1. Lange, Isabelle, and S. E. Forbush, Cosmic-ray results 
from Huancayo Observatory, Peru, June 1936- 
December 1946, Carnegie Inst. Wash. Publ. 175, 
vol. XIV, Washington, D. C., 1948. 

2. Forbush, S. E., Time-variations of cosmic rays, Hand- 
buch der Physik. Band XLIX/I, Geophysik III, Teil 
I, pp. 159-247, Springer-Verlag, Berlin, Heidelberg, 
New York, 1966. 

3. Forbush, S. E., and D. Venkatesan, J. Geophys. Res., 
Vol. 65, No. 8, August 1960, pp. 2213-2226. 


DESCRIPTION OF TABLES 


INSTITUTO GEOFISICO DE HUANCAYO 


Tables 1-108. Monthly tables for the period January 
1960 through December 1968 contain bihourly mean 
values of departures from balance corrected to 510 mm 
Hg. Bihourly mean values are centered at the indicated 
Greenwich mean hours. Tabular values are in units of 0.1 
percent of total cosmic-ray ionization with 6.0 percent 
added. 

Tables 109-117. Annual tables for the period 1960 
through 1968 contain daily means of hourly departures 
from balance, corrected to 510 mm Hg. Tabular values, 
by months, are in units 0.1 percent of total cosmic-ray 
ionization with 6.0 percent added to obtain the tabular 
values. Superscripts a to k indicate the number of missing 
bihourly values | to I1. 

Tables 118-126. Annual tables for the period 1960 
through 1968 contain bihourly values of diurnal vari- 
ation in departures from balance, corrected to 510 mm Hg, 
centered at the indicated Greenwich mean hours. Tabu- 
lations are by months, seasons, and years, from days with 
complete records. Means of three groups of four months 
each are designated by roman numerals as follows: I, 
January, February, November, December; II, March, 
April, September, October; III, May, June, July, August. 
Tabular values are in units of 0.01 percent of total cosmic- 
ray ionization, with 6.0 percent added. 


FREDERICKSBURG MAGNETIC OBSERVATORY 


Tables 127-234. Monthly tables for the period Jan- 
uary 1960 through December 1968 contain bihourly mean 
values of departures from balance, corrected to 760 mm 
Hg. Bihourly mean values are centered at the indicated 
Greenwich mean hours. Tabular values are in units of 0.1 
percent of total cosmic-ray ionization, with 10.0 percent 
added. 

Tables 235-243. Annual tables for the period 1960 
through 1968 contain daily means of hourly departures 
from balance, corrected to 760 mm Hg. Tabular values, 
by months, are in units of 0.1 percent of total cosmic-ray 
ionization with 10.0 percent added to obtain the tabulated 
values. Superscripts a to k indicate the number of missing 
bihourly values | to 11. 

Tables 244-252. Annual tables for the period 1960 
through 1968 contain bihourly values of diurnal variation 
in departures from balance. corrected to 760 mm Hg, 
centered at the indicated Greenwich mean hours, by 
months, seasons, and years, from days with complete rec- 
ords. Means for three groups of four months each are 


designated by roman numerals as follows: I, January, 
February, November, December; II, March, April, 
September, October; HI, May, June, July, August. 


Tabular values are in units of 0.01 percent of total cosmic- 
ray ionization, with 10.0 percent added to obtain the tabu- 
lated values. 


CHRISTCHURCH MAGNETIC OBSERVATORY 


Tables 253-255. Annual tables for the period 1959 
through 1961 contain daily means of hourly departures 
from balance, corrected to 29.00 in. Hg. Tabular values, 
by months, are in units of 0.1 percent of total cosmic-ray 
ionization. Superscripts a to k indicate the number of 
missing bihourly values 1 to 11. 

Tables 256-258. Annual tables for the period 1959 
through 1961 contain bihourly values of diurnal variation 
in departures from balance, corrected to 29.00 in. Hg, 
centered at the indicated Greenwich mean hours, by 
months, seasons, and years, from days with complete rec- 
ords. Means for three groups of four months each are 
designated by roman numerals, as follows: I, January, 
February, November, December; I, March, April, 
September, October; III, May, June, July, August. Tabular 
values are in units of 0.01 percent of total cosmic-ray 
ionization. 


GODHAVN MAGNETIC OBSERVATORY 


Tables 259-294. Monthly tables for the period Jan- 
uary 1951 through December 1953 contain bihourly mean 
values of departures from balance corrected to 760 mm 
Hg. Bihourly mean values are centered at the indicated 
Greenwich mean hours. Tabular values are in units of 0.1 
percent of total cosmic-ray ionization, with 13.0 percent 
added. 

Tables 295-297. Annual tables for the period 1951 
through 1953 contain daily means of hourly departures 
from balance, corrected to 760 mm Hg. Tabular values, 
by months, are in units of 0.1 percent of total cosmic-ray 
ionization, with 13.0 percent added to obtain the tabulated 
values. Superscripts a to k indicate the number of missing 
bihourly values | to 11. 

Tables 298-300. Annual tables for the period 1951 
through 1953 contain bihourly values of diurnal variation 
in departures from balance, corrected to 760 mm _ Hg, 
centered at the indicated Greenwich mean hours, by 
months, seasons, and years, from days with complete rec- 
ords. Means of three groups of four months each are 
designated by roman numerals, as follows: I, January, 
February, November, December; II, March, April, 
September, October; HI, May, June, July, August. 
Tabular values are in units of 0.01 percent of total 
cosmic-ray ionization, with 13.0 percent added to obtain 
the tabulated values. 
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